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ABSTRACT 

We search for an unusual alignment of the preferred axes of the quadrupole and octopole, 
the so-called axis of evil, in the CMB temperature and polarization data from WMAR We use 
the part of the polarization map which is uncorrelated with the temperature map as a statis- 
tically independent probe of the axis of evil, which helps to assess whether the latter has a 
cosmological origin or if is a mere chance fluctuation in the temperature. Note, though, that 
for certain models creating a preferred axis in the temperature map, we would not expect to 
see the axis in the uncorrelated polarization map. We find that the axis of the quadrupole of 
the uncorrelated polarization map roughly aligns with the axis of evil within our measure- 
ment precision, whereas the axis of the octopole does not. However, with our measurement 
uncertainty, the probability of such a scenario to happen by chance in an isotropic universe 
is of the order of 50 per cent. We also find that the so-called cold spot present in the CMB 
temperature map is even colder in the part of the temperature map which is uncorrelated with 
the polarization, although there is still a large uncertainty in the latter. Therefore, our analysis 
of the axis of evil and a future analysis of the cold spot in the uncorrelated temperature data 
will strongly benefit from the polarization data expected from the Planck satellite. 
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1 INTRODUCTION 

A major assumption of modern day cosmology is the cosmologi- 
cal principle, which states that the Universe is homogeneous and 
isotropic on large scales. The observed isotropy of the Cosmic Mi- 
crowave Background (CMB) is one of the strongest evidences sup- 
porting the cosmological principle. 

However, in recent years, there have been claims of anoma- 
lies detected in the CMB temperature map with consider- 
able significance, which seem to break statistical isotropy of 
the temperature fluctuations and thus to question the cosmo- 



ferred direction in the CMB temperature map have led to a dis- 
cussion about whether this is simply due to a chance fluctuation 
in the CMB temperature map, if it can be blamed on local struc- 
tures or on systematics in the measurement, or whether it is actually 
due t o a preferred direct i on intrinsic to our Universe jCopi et al.1 
| 200l iDolag et al.1 120051: iMaturi et al.1 [20071: ISamal et all l2009r 
iGroeneboom & Erik sen 2009: iMorales & Saez 2008; Vielva et al 
20071: llnoue& Silkl [20071: [Gaol [20091: lAckerman et al.1 | 2007| 
Cooi et al.1 120061: ISchwarz et al.1 120041 : lHansen et al.1 12004 2009 



Prunet et al.1 120051 : Ijaffe et al.1 1200I 120061 : iBernui et all l20O 
Wiaux et al.ll2006h . 



logical principle. Several groups ( de Oliveira-Costa et al 
Abra mo et alJ l200d: lLand & Magueiiol l200l I Samal etal 



2004; 



2008 



Raki c & Schwarzl 120071) claim to have found a strong alignment 



between the preferred axes of the quadrupole and the o ctopole, 
which is commonly referred to as the axis of evil. Others dBernuil 
l2008l : lEriksen et alfeOOllHoftuft et alJl2009h have found a signif- 
icant power asymmetry between the northern and southern eclip- 
tic hemisphere, and some weaker anomalies have been found 
for the low multipoles beyond the octopole dCopietalJ 12 004; 
Land & Magueiiol 120051 : lAbramo et all 120061 : IPereira & Abramol 
2009). However, the existence of such an isotropy breaking 



in the CMB temperature map is strongly under debate, and 
also negative results have been published dSouradeep et al.l 1 2006; 
iMagueiio & So rkin 2007). The claims of the existence of a pre- 



The polarization fluctuations of the CMB, just as its tempera- 
ture fluctuations, have their origin in the primordial gravitational 
potential. The polarization should thus exhibit similar peculiari- 
ties as the temperature, provided they are due to some preferred 
direction intrinsic to the geometry of the primordial Universe. Note 
that this is not generic to every theoretical model creating anoma- 
lies in the temperature map. For example, if the peculiarities in 
the temperature maps are due to a secondary effect on the CMB 
such as the integrated Sachs- Wolfe effect, we would not expec t 
them to be present in the polarization maps tovorkin et al.l l2008). 
The search for anomalies in the CMB polarization map is still in 
its initial stage, due to the high noise-level in the available full- 
sky polarization map from the W ilkinson Microwave Anisotropy 
Probe (WMAP). ISouradeep et al.l ( 2006) have found some evi- 
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Table 1. Axes and their uncertainties for the four different maps in Galactic 
coordinates. The large errors are due to the effects of the mask, residual 
foregrounds and the detector noise in the WMAP polarization data. 



dence for anisotropi es in the WMAP p olarization data using the 
method proposed in Bas ak et all d2006h . However, they state that 
the anisotropies are likely due to observational artifacts such as 
foreground residuals, and that further work is required in order to 
confirm a possible cosmic origin. 

Given that the polarization map is correlated with the temper- 
ature map, it is not a statistically independent probe of the anoma- 
lies which have been found in the temperature map. If the ob- 
served anomalies were due to a chance fluctuation in the tempera- 
ture map, this chance fluctuation could also be prese nt in the polar- 
izatio n maps, due to the correlation between the two dDvorkin et ajj 
2008). In this work, we split the WMAP polarization map into a 
part correlated with the temperature map, P CO rr, and a part uncor- 
rected with the latter, Puncorr- We obtain the part of the polariza- 
tion map which is correlated with the temperature map by simply 
translating the temperature map into a polarization map, using their 
cross-correlation. The part of the polarization map which is uncor- 
rected with the temperature map serves as a statistically indepen- 
dent probe of the above-mentioned anomalies. Chance fluctuations 
in the temperature maps do not affect the uncorrected polarization 
map, so that a detection of the anomalies in the latter would be a 
hint to an actual cosmological origin of them. Note, though, that 
this does not have the power to exclude residual foregrounds or 
systematics as potential origins for the anomalies. 

Similarly, we split the WMAP temperature map into a part 
correlated with the polarization map, T^r, and an uncorrelated 
map, Tuncorr- If the anomalies detected in the CMB temperature 
map are of genuine cosmological origin, they should be present in 
the correlated and the uncorrelated parts of both the temperature 
and polarization map. For convenience, the four resulting maps are 
summarised and briefly described in Table Q] 

In this work, we focus on using the uncorrelated polariza- 
tion map to probe the axis of evil. In order to define the pre- 
ferred axis of the multipole s, we use a statistic proposed by 
Ide Oliveira-Costa et al.l ([2004), which is the axis around which the 
angular momentum dispersion is maximised for a given multipole /. 
We note that we have to mask out about 25 per cent of the sky in the 
WMAP polarization data in order to reduce Galactic foregrounds. 
Furthermore, the polarization data are highly contaminated by de- 
tector noise and residual foregrounds even outside the mask. We 
therefore perform a Wiener filtering of the polarization data before 
determining the preferred axes, in order to reduce the noise con- 
tained in the maps. However, we still expect a large uncertainty in 
the axes, which we obtain by running Monte Carlo (MC) simula- 
tions conditional on the data. The uncertainty in our axes amounts 
to cr^45°. 

We find that, for all four of the maps, the preferred axes of 
the quadrupole all point in the same direction, within our measure- 



ment precision. However, the preferred axis of the octopole of the 
uncorrelated polarization map does not align with the one of the 
quadrupole. The same holds for the correlated temperature map. 

In order to assess our result, we ask the following question. 
We take the axes measured in the temperature map as given, and 
assume that the axes of the uncorrelated polarization map are dis- 
tributed isotropically and independently of each other. We then ask 
how likely it is that at least one of these axes lies such that the axis 
of the temperature map lies inside its la region. This probability 
amounts to about 50 per cent for currently available polarization 
data. This high probability is due to the large uncertainties we have 
in the axes of the uncorrelated polarization map. The main contri- 
bution of this uncertainty comes from the high noise-level in the 
polarization data rather than from t he mask. We c an therefore hope 
that the Planck polarization data (Tauberi l2000h will yield much 
stronger constraints on the axes than the WMAP data. 

Note that our approach to probing the axis of evil in polar- 
ization is phenomenological, since not all theoretical models of 
the primordial Universe exhibiting anomalies in the CMB tempera- 
ture map show the same behaviour in the uncorrelated polarization 
map. We outline a more thorough analysis, taking into account the 
predictions of the specific models for the uncorrelated polarization 
map, in the conclusions of this work. 

The article is organised as follows. In section [2 we briefly re- 
view the Wiener filter. In sections [3] and HI we explain in detail the 
splitting of the WMAP temperature and polarization maps, respec- 
tively. Section[5]is devoted to determining the preferred axes for the 
quadrupole and octopole for our four maps. We conclude in section 

m 



2 WIENER FILTERING 

As we have already mentioned in the introduction, the WMAP 
polarization data are highly contaminated by detector noise and 
Galactic foregrounds. The observed polarization map we use is the 
linear combination of the maps of the Ka, Q, and V frequency bands 
(corresponding to 33, 41, and 61 GHz), which is used for determin- 
ing the lo w-1 polarization likelih ood in the 5 year WMAP likeli- 
hood code jHinshaw et al.ll2009b . By using the linear combination 
of the maps, we combine the information from different frequency 
bands rather than using only the information contained in a particu- 
lar band. Therefore, the linear combination is less contaminated by 
noise than the original maps per frequency band. We use the P06 
mask JPage et aDl2007h to mask out the Galactic plane in the polar- 
ization map. The linear combination maps for the Stokes Q and U 
parameters are shown in Fig.Q]in Galactic coordinates. 

In order to reduce the noise level, we perform a Wiener filter- 
ing of the observed polarization map before translating it into the 
part of the temperature map which is correlated with the polariza- 
tion data. Similarly, we will perform a Wiener filtering of the part 
of the polarization map which is uncorrelated with the temperature 
map, as we will describe in detail later on. 

The Wiener filter can be derived for the general data model 

d = FLs + n, (1) 

where d denotes the data, s the (temperature or polarization) signal, 
R the instrument response and n additive noise. Let us define the 
signal and noise covariances, 

S = (ss f ) P(s) = [ Vs{ss ] )V(s), (2) 
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Figure 1. Observed polarization maps (linear combination of Ka, Q, and V 
band maps). Stokes Q map (top panel) and Stokes U map (bottom panel). 
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where the dagger denotes a transposed and complex conjugated 
quantity, V(s) and V(ri) denote the probability density functions 
of signal and noise, respectively, and the integrals have to be taken 
over all pixels i, e.g. 



(4) 



If we assume the signal prior and the noise distribution to be Gaus- 
sian, we obtain the signal posterior 



V(s\d) = g(s-s Tec ,D) 

Here, we have defined 
1 



Q(x,C) 



exp 



(5) 



(6) 



to denote the probability density function of a Gaussian distributed 
vector x with zero mean, given the cosmological parameters p and 
the covariance matrix C = (xx*)» where the averages are taken 
over the Gaussian distribution G(x, C). In eq. (|5]), we have used 
the definitions 

srec = (S' 1 + R i N~ 1 R)~ 1 R i N~ 1 d , (7) 
which is called the Wiener reconstruction of the signal, and 

D = (S' 1 + R i N~ 1 R)~ 1 , (8) 

which denotes the Wiener variance with which the real signal, s, 
fluctuates around the reconstruction, s re c. A detailed derivation 
of the posterior dis tribut ion, eq. ([5]), can for ex ample be found in 
lEnfilin et al.1 (120081) or in lFrommert etaD <2008h . 



3 SPLITTING OF THE TEMPERATURE MAP 

In this section, we split the WMAP temperature map into a part cor- 
related with the WMAP polarization map, T^r, and a part which 
is not, T*n C coTT- This is the same splitting which has been done in 
iFrommert & EnBlinl (12009 ) in order to reduce the noise in ISW 
measurements. We translate the polarization map into the corre- 
lated part of the temperature map, using the cross-correlation be- 
tween the two. However, as we already mentioned in the last sec- 
tion, before doing so we perform a Wiener filtering of the observed 
polarization map in order to reduce the noise. 

Our data model for the observed polarization map P, which 
contains the Stokes Q and U maps shown in Fig.Q] is 



Q 
u 



W (P cmb + Pdet + Pf s 



(9) 



Here, P cm b is the intrinsic CMB polarization, Pdet and Pf g denote 
the detector noise and residual foregrounds, respectively, and we 
have introduced the window W in order to describe the mask. 

Let us define the signal covariance matrix of the CMB polar- 
ization given the cosmological parameters p, 



(Per 



t\ 



(10) 



and the noise covariance matrices for the detector noise and the 
residual foregrounds: 



N det 
JVfir 



t) > 



(Pdet Pdet ' )V(P de 

(Pf gJ Pfg t )p(P fg ) . (ID 

The signal power spectru m (and thus Sp) has been computed using 
CMBEASY fooranll2005l) for the Maxi mum Likelihood cosmolog- 
ical model from lDunklevetal.1 d2009h : {Q b h 2 = 0.0227, Q A = 
0.751, h = 0.724, r = 0.089, n s = 0.961, a 8 = 0.787}. 

In order to derive the Wiener filter for P, let us define the 
noise, 



Tl = W (Pdet + Pfg) , 

for which the noise covariance is then 

N P = {nn ] ) P(n) = W(N det + N (s )W f , 



(12) 



(13) 



where we have assumed that Pdet and Pf g are uncorrected. We take 
the total noise covariance, Np, for the observed polarization map 
from the WMAP code. We further identify P cm b with the signal 
s, the mask W with the response R, and P with the data d. With 
these definitions, we have translated our data model, eq. (|9]), into 
the one given in eq. ([TJ. If we assume the noise n and the signal 
Pcmb to be Gaussian distributee^ we therefore obtain the posterior 
distribution for the signal 

P(Pcmb | P,P) = G (Pcmb - Pcmb, AO , ( 14 ) 

with 

Pcmb = (Sp 1 + W^N^Wy^N^P , (15) 
which is the Wiener reconstruction of the polarization map, and 

D p = (Sp 1 + W^N^Wy 1 , (16) 

which denotes the Wiener variance. We show the Stokes Q and 
U maps of the Wiener filtered polarization map Pcmb m me t0 P 

1 The assumption of Gaussianity holds well for the detector noise Pd e t and 
the signal P cm b- F° r the residual Galactic foregrounds, this assumption is 
probably less accurate. 
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Figure 2. Stokes Q part of the following polarization maps: Top panel: 
Wiener filtered polarization map, P?^ b . Middle panel: Part of the polariza- 
tion map correlated with the temperature map, P CO rr- Bottom panel: Part 
of the polarization map uncorrelated with the temperature map, Puncorr- 
The colour scale is the same in all maps. 



panels of Fig. |2] and Fig. [3] respectively. Note that only the low / 
modes survive the Wiener filtering, whereas the higher / modes are 
strongly suppressed due to the high noise-level they contain. 

We now split the WMAP temperature map into a part cor- 
related with the polarization map, Tco? r » and a part uncorrelated 
with the latter, Puncorr- We use the Wiener filtered polarization 
map Pcmb» which is of resolution NSID E=8, and the inte rnal lin- 
ear combination (ILC) temperature map ( Gol d et al . 2009), which 
we have smoothed with a Gaussian beam of FWHM=18.3° and 
downgraded to the same resolution. Among the different WMAP 
temperature maps, the ILC is the one for which the alignment 
of the lo w multipoles is least conta minated by Galactic fore- 
grounds dGruppuso & Buriganal 120091) . When working on large 
scale s, we can safely neglect the detector noise in the temperature 
data (Afshordi 2004). Furthermore, we decide to neglect residual 
foregrounds in the temperature map. 

We translate the Wiener filtered polarization map, P c r m b , 
into the correlated part of the temperature map, using the cross- 
correlation between the two: 

T?rec o q— 1 D rec f\n\ 
corr = &T,P Op r C mb j U ') 




Figure 3. Stokes U part of the following polarization maps: Top panel: 
Wiener filtered polarization map, P5^ b . Middle panel: Part of the polariza- 
tion map correlated with the temperature map, P CO rr- Bottom panel: Part 
of the polarization map uncorrelated with the temperature map, -Puncorr- 
The colour scale is the same in all maps. 

where the signal covariance matrices given the cosmological pa- 
rameters, p, are defined as 

S P ,T = (PcmbT t )p (T)Pcmb | p) , (18) 

St = (TT^) viT]p) . (19) 

The uncorrelated temperature map T^corr is then obtained by 
simply subtracting T^r from T: 

-L uncorr — -*■ -*■ corr • 

In Appendix lAl we prove that and T^corr are indeed uncor- 
related. 

We plot T, T c r orr, and T^corr m the top, middle, and bottom 
panel of Fig.|4] respectively. Let us first concentrate on T^r, and 
try to assess whether some of its structures could come from Galac- 
tic foregrounds rather than being intrinsic CMB fluctuations. Note 
that this is just meant to be a quick glance on what we can imme- 
diately pick out by eye. Comparing P^orr with the overview over 
the Galactic foregrounds published in lHinshaw et al.l d2007h . Fig. 
7, makes us suspect that the warm region in the middle of the north- 
ern hemisphere might be associated with the North Galactic Spur. 
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Figure 5. The cold spot, which has been found to have non-Gaussian char- 
acteristics, is marked in the ILC map shown here by a black circle. 





Figure 4. Top panel: ILC map, smoothed with a beam of 18.3° and down- 
graded to a resolution of NSIDE=8. Middle panel: Part of the temperature 
map which is correlated with the polarization map, T™£ T . Bottom panel: 
Part of the temperature map which is uncorrected with the polarization 
map, T^ncorr- The colour scale is the same in all maps. 



A part of this region is already masked out, but it is well possi- 
ble that the mask should be bigger in order to better mask out this 
foreground. One might also think that the big red blob on the right 
hand side of T^r, close to the Galactic plane, could be due to the 
Gum Nebula. However, plotting the two maps on top of each other 
reveals that the Gum Nebula lies further to the East than our red 
blob. Therefore we exclude that the blob comes from that particu- 
lar foreground. 

Let us now compare the three maps T, T^r, and T^corr- m 
the northern Galactic hemisphere, all maps look quite similar, apart 
from the hot region around the North Galactic Spur, which is more 
prominent in T^rr than in the other two maps, and which we have 
already commented on. However, in the western part of the south- 
ern hemisphere, we obtain a strong deviation of from the ILC 
map. In fact, the features in T££ T have the opposite sign to the 
structures in the ILC map. This enhances the amplitudes of the fea- 
tures in the western part of the southern hemisphere in T^corr as 
compared to the ILC map. In particular, the so-called cold spot, 
which has been fo und to have non-Gaussian characteristics by 
IVielva et all (120041) , turns out to be even colder in T^corr than in 
the ILC map. The cold spot, which we mark in the ILC map in Fig. 



\5\ by a black circle, has later been confir med to have non-Gaussian 
characteristics by many others (see, e.g.jMartmez- Gonzalez et all 
2006: ICruz eta l. 2006; Naselsky et al. 200^ It would be interest- 
ing to redo the above-mentioned analyses of the cold spot with the 
high-resolution version of Tuncorr> m order to see whether the sig- 
nificance of the non-Gaussian features is even higher in that map. 
A thorough analysis of the characteristics of the cold spot is beyond 
the scope of this work, though, and we leave this exciting question 
for future work. Lastly, we notice that on the large scales we are 
looking at, we have much stronger deviations of the temperature 
towards the cold end of the temperature spectrum than towards the 
warm end, for all three of the maps. 



4 SPLITTING OF THE POLARIZATION MAP 

We now split the WMAP polarization map into a part correlated 
with the WMAP temperature map, P CO rr, and a part uncorrected 
with that, Puncorr- As before, we obtain the correlated polarization 
map by simply translating the temperature map into a polarization 
map: 



Pcorr = SpsrSrp^T , 



(21) 



The Stokes Q and U maps of P CO rr are shown in the middle panels 
of Fig.|2]and Fig. [3] respectively. 

In order to obtain the uncorrelated map, we would like to sub- 
tract Pcorr from P cmb : 



= P C1 



(22) 



However, we do not know P cm b because we only observe P, which 
is highly contaminated by noise. Subtracting P CO rr from the Wiener 
filtered polarization map, Pcmb» does not result in uncorrelated 
maps. We therefore subtract WP CO tt from the observed polariza- 
tion map, P: 



P - W Pcorr 
W Puncorr + 71 , 



(23) 



where the noise n is the same as in section [3] We then compute 
the Wiener reconstruction of the signal Puncorr, with the data being 

nraw 
uncorr • 

Puncorr = [{S P - S P ,tS^ St.p) + VK f Np V] 

W^Np'P^corr- (24) 

Here, we have used the signal covariance 

(Puncorr Puncorr ^)p(P cmb ,T | p) 
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— (-Pcmb-Pcmb^) — (PcmbP^ ) S T 1 St,P 

— Sp^rSrp 1 (TP cm b^) + Sp^rSrp 1 (TT^)Srp 1 St,p 

— Sp — Sp^tSji 1 St,p • (25) 

Puncorr given in eq. <f24t is uncorrected with P CO rr, as we prove in 
Appendix iBl The posterior of Puncorr is given by 

ncorr I P 5 P, p) — Q (Puncorr Puncorr 5 A, ncorr ) , (26) 

with the Wiener variance 

Amcorr = [(S P ~ S P ,T 5 T ,p) _1 + W 1 N^W]' 1 . (27) 

We show the Stokes Q and U maps of the uncorrected po- 
larization map, Puncorr* m the bottom panels of Fig. [2] and Fig. [3] 
respectively. Note that the symbols for the correlated and uncor- 
rected parts of temperature and polarization maps are listed and 
briefly explained in Table Q] 



5 THE AXIS OF EVIL 

We now search for the axis of evil in the four maps P CO rr, Puncorr. 
T c r orr, and T^corr- Note that P C orr and T^rr have of course the 
same axes as the original temperature and polarization maps, T 
and Pcmb, resp ectively. To define the preferred axis, we use a statis- 
tic proposed bv lde Oliveira-Costa et al.1 12004), which has been in- 
troduced in order to quantify the preferred direction that can be 
picked out in the smoothed temperature map by eye. When look- 
ing at the smoothed ILC map in Fig. |H most of the hot and cold 
blobs seem to be lying on the same plane. The quadrupole and oc- 
topol e extracted from the ILC map show the same behaviour (see, 
e.g.. Ide Oliveira-Costa et al.|[2004l) . and the planes are roughly the 
same for the two multipoles. In order to quantify this alignment, 
Ide Oliveira-Costa et"aD J2004) came up with the following statis- 
tic. The temperature maps are expanded into spherical harmonics, 
which are eigenfunctions of the square and the z-component of the 
angular momentum operator L: 

T(n) = ^T z (n) = ^aL^(n) . (28) 

l l,m 

Then, for every multipole /, one determines the z-axis n for which 
the expectation value of the z-component of L, ft L, is maximised: 

(T, | (ft • Lf | T,) = J2 ™ I a Tm(n) | 2 , (29) 

m 

Here, af m (n) denotes the spherical harmonic coefficient af m ob- 
tained in a coordinate system with the z-axis pointing in re- 
direction. We determine the axis ft by simply rotating the z-axis 
into every pixel centre and checking for the maximum, which is 
well feasible at our resolution. Neighbouring pixel centres in our 
map differ by approximately 7° , but we will soon see that the un- 
certainties in our axes are so large that it is sufficient to check only 
the pixel centres as potential z-axes. We have done the same ex- 
ercise allowing the axes to point to all pixel centres of NSIDE=16 
instead of NSIDE=8, and our results are robust under this change. 

As we have already mentioned, the mask, residual foregrounds 
and detector noise in the polarization data will result in an uncer- 
tainty in the preferred axes. The posterior distribution of the real 
CMB polarization map, P cm b, given the one we observe, P, is 
given by eq. dT4t . P cm b fluctuates around our Wiener reconstruc- 
tion, Pcmb» w i m the Wiener variance Dp. 




Figure 6. Wiener realisation of T cor r 

In order to obtain the uncertainties in the axes of T™* T and 
Tu n corr> we have run Monte Carlo (MC) simulations, drawing re- 
alisations of Pcmb from its posterior distribution. From these, we 
obtain realisations of 

T^corr = St,P Sp 1 Pcmb •> 
Puncorr = T Pcorr •> (30) 

for which we then determine the preferred axes. The uncertainty in 
the axes of Puncorr is obtained similarly, using the posterior dis- 
tribution of Puncorr given in eq. d26b . Note that T and thus Pcorr 
are assumed to have no contributions from residual foregrounds or 
detector noise, and thus no uncertainty in the preferred axes. 

For drawing realisations from the probability distribution in 
eq. dl4l , we have computed the Wiener variance Dp given in eq. 
JTol . We have then computed the Cholesky factorisation L of Dp, 
which is a particular form of the square-root of a positive definite 
matrix: 

D P = LL t . (31) 

In order to obtain our realisation, P cm b, we apply L to a map n w of 
white noise, i.e. a map where the temperature at every pixel is inde- 
pendently drawn from a Gaussian distribution with unit variance, 
and add the mean value Pcmb : Pcmh = Ln w + Pcmb- This results 
in a map which is drawn from the distribution in eq. (fl4l . as one 
can easily see: 

((Pcmb — Pcmb) (-Pcmb — P C mb ) ^ ) V(n w ) 

- L (n w nl,)<p( nw ) P f = PP f = Dp. (32) 

An example of a Wiener realisation of T CO rr in shown in Fig. [60 

We plot the axes and their uncertainties for the different maps 
in Figs|7]-[T0] Both ends of every axis are marked by a cross in the 
maps, and the colour coding counts how many times the preferred 
axis came to lie on the respective pixels in 5000 MC samples. 

All axes and their standard deviations a, which we ob- 
tained from the MC simulations, are summarised in Table [T] 
For Pcorr, and thus the ILC m ap, we reproduce the results from 
Ide Oliveira-Costa et al.l |2004|) within our measurement precision: 
the axes of the quadrupole and the octopole of P CO rr point in 
the same direction, which is roughly (Z,6) « ( — 120°, 63°), 

2 We had to regularise the Wiener variances, eqs ( fl6l and {27), by adding 
Gaussian noise in order to make them positive definite. This is required by 
the Cholesky factorisation. However, since the noise was added mostly on 
small scales, the quadrupole and octopole remained completely unaffected 
by this. In fact, our results remained unchanged under varying the variance 
of the added Gaussian noise over 5 orders of magnitude. 
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Figure 7. Preferred axis of the quadrupole (top panel) and the octopole 
(bottom panel) for P CO rr and thus for the ILC map. We reproduce the results 
of Ide Oliveira-Costa et al.l d2Q04l) within our measurement precision. The 
axes of quadrupole and octopole point in the same direction, which has 
been named the axis of evil. 

where / and b denote Galactic longitude and latitude, respectively 
(Ide Oliveira-Costa et all (120041) found (Z,6) « (-110°, 60°)). For 
Tuncorr> again both axes point in the same direction as the axes of 
Pcott within our measurement precision. 

For Puncorr, the preferred axis of the quadrupole has an angu- 
lar distance to the average axis of the ILC map of 37° . That means 
that the latter lies inside its la region. The same holds for 
(and thus PcmbX for which the axis of the quadrupole has an angu- 
lar distance to the average axis of the ILC map of 34° . The axes of 
the octopole of Puncorr and T^r, though, do not align with the axis 
of evil. What can we learn from this result? The significance of the 
alignment between the axes of the quadrupole and octopole in the 
temperature map has been assessed extensively in earlier works. 
In this work, we only look at the additional information we ob- 
tain from the axes of Puncorr- To this end, let us take the preferred 
axis in the temperature map T as given, and assume that the axes of 
Puncorr are distributed isotropically over the sky and independently 
from each other. In Appendix [Cl we work out the probability for at 
least one of the axes of Puncorr being such that the axis of the tem- 
perature map is included in the la region around it. This probability 
amounts to about 50 per cent, due to the large la regions we have. 

In order to assess whether the mask or the noise in the WMAP 
polarization maps is the main source of uncertainty in the axes, we 
have determined the uncertainty with the amplitude of the noise 
covariance matrix rescaled to 10 per cent of the original one. This 
yields an uncertainty of about 20° in the axes. We have done the 
same exercise for the noise amplitude downscaled to 1 per cent of 
the original one, which results in an uncertainty of 7° — 10° in 
the axes. This means that the noise is actually the main source of 
uncertainty in our analysis rath er than the mask. Soon, the Planck 
Surveyor mission (TaubeJ bOOOh will provide us with polarization 
measurements that have a noise-level which is significantly below 
the one in the WMAP data. The main problem will then be the 
contamination of the polarization data by Galactic foregrounds. In 
the WMAP polarization data, the foregrounds contribute about 20 
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Figure 8. Preferred axis of the quadrupole (top panel) and the octopole 
(bottom panel) for Puncorr- The colour coding counts the number of MC 
samples whose axis came to lie on the respective pixel. The axis of the 
quadrupole aligns with the axis of evil within our measurement precision, 
whereas the axis of the octopole does not. 




Figure 9. Preferred axis of the quadrupole (top panel) and the octopole 
(bottom panel) for T^ r and thus for P^ h - The axis of the quadrupole 
aligns with the axis of evil within our measurement precision, whereas the 
axis of the octopole does not. 
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Figure 10. Preferred axis of the quadrupole (top panel) and the octopole 
(bottom panel) for T^ COTT . The axes of the quadrupole and the octopole 
both align with the axis of evil within our measurement precision. 

per cent to the diagonal of the noise covariance matrix Np in pixel 
space. With Planck, we will be able to determine the foregrounds 
better than with WMAP, due to the broader frequency range cov- 
ered by Planck. If we assume that the covariance due to residual 
foregrounds for Planck will be between 5 and 50 per cent of the 
one for WMAP, we will get the uncertainty on the axes down to 
about 10° and 20°, respectively. With this, we will have a powerful 
test to probe the axis of evil in polarization. 



6 CONCLUSIONS 

In the last few years, a preferred axis has been found in the CMB 
temperature map, posing a challenge to the cosmological principle. 
This so-called axis of evil denotes the unusual alignment of the pre- 
ferred axes of the quadrupole and the octopole in the temperature 
map. 

In this work, we have split the CMB temperature and polar- 
ization maps from WMAP into a part correlated with the respective 
other map, and an uncorrelated part. If the axis of evil were due to 
some preferred direction intrinsic to the geometry of the primordial 
Universe, we would expect its signature to be present in all four of 
these maps, although this is not true for all theoretical models cre- 
ating an axis in the temperature map. In particular, the part of the 
polarization map which is uncorrelated with the temperature map 
serves as a statistically independent probe of the axis of evil. 

In order to reduce the noise contained in the polarization maps, 
we have Wiener filtered the maps before computing the axes. We 
have then determined the preferred axes of the quadrupole and the 
octopole for the four maps. In order to assess the uncertainty in 
the axes coming from the mask, detector noise and residual fore- 
grounds in the polarization maps, we have run MC simulations con- 
ditional on the observational data. 

For the part of the polarization map which is correlated with 



the temperature map, P CO rr, we find that the axes of quadrupole 
an d octopole point in the same direction, confirming earlier results 
bv lde Oliveira-Costa et al.1 j2004l) . The part of the temperature map 
which is uncorrelated with the polarization map, T^ corr , exhibits 
the same alignment of the axes within our measurement precision. 

For the part of the polarization map which is uncorrelated 
with the temperature map, Puncorr> we find that only the axis of 
the quadrupole aligns with the axis of evil, whereas the axis of the 
octopole does not. The same holds for the correlated part of the 
temperature map, T^r- We have computed the probability that a 
rough alignment with the axis of evil, as we find it for the axis of 
the quadrupole of Puncorr* happens by chance if the axes are dis- 
tributed isotropically. This probability amounts to 50 per cent for 
currently available polarization data, due to the large uncertainties 
in the axes. We are thus looking forward to redoing this analysis 
with polarization maps from Planck, which will yield much more 
significant results. Of course, similar analyses can be carried out for 
all other anomalies that have been found in the CMB temperature 
maps. Note that, instead of working in pixel space as we have done, 
one could implement the analysis in spherical harmonics space, 
which would help to separate the E modes we are working with 
from contamination by B modes. 

The approach we have chosen here is a phenomenological ap- 
proach, since in principle one should take into account that different 
models causing anomalies in the temperature map predict differ- 
ent signatures in the polarization map. Thus, for a more thorough 
analysis, one would need to consider particular models of the pri- 
mordial Universe creating anomalies in the temperature maps, and 
compute the statistical properties of the uncorrelated polarization 
map for these. This can be done by modifyi ng a Boltzma n n cod e 
such as CMB EASY or by simulations as in iDvorkin et al.1 d2008h . 
One can then try to find these predicted signatures in the uncorre- 
lated polarization map via Bayesian model selection. Such an anal- 
ysis would truly go beyond the usual a posteriori analysis of finding 
anomalies in the temperature map, since we would use an actual 
model to make predictions for the uncorrelated polarization map 
and then compare these predictions with observations. We leave 
this promising analysis for future work. 
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APPENDIX A: PROOF OF VANISHING CORRELATION 
BETWEEN T c r e r c r AND T u T corr 

We now prove that the two maps and T^ c corr , into which we 
split the temperature map T in section [3] are indeed uncorrected. 
To this end, let us write 

corr — £>T,P Op -*cmb 

= Sr^Sp^Sp 1 + Np l W)~ l W 1[ Np 1 P 

S t ,p(1 + W^N^WSpy^^N^P 
= S T ,pW\l + N^WSpW^^N^P 
= S^pW^Np + WSpW^P, (Al) 

where we have inserted P c r mb from eq. (1131) in the first step. The 
third step can be easily verified by using the geometric series for 
(l + W^N^WSp)' W\ which has a convergence radius of 1, 
and is thus valid for \W^ N~ 1 WS P \ < 1. In our case, this holds 
because our polarization data are noise-dominated[f| 

We will soon see that we need the following covariance matri- 
ces in the derivation: 

(PP t ) P(P | p) = WS P W ] + Np , (A2) 

where we have assumed that P cm b is uncorrelated with Pdet and 
Pf g , and we have inserted the definition of Np, eq. ( fT3l ). 

Since we neglect the detector noise and residual foregrounds 
in the temperature data, we obtain for the covariance between tem- 
perature and polarization data 

(TP*)vtT,P\ P ) = (TPcmb t )p(T J P cmb |p)W t 

= S T ,pW* , (A3) 

where we have assumed that detector noise and residual fore- 
grounds in the polarization map are uncorrelated with the CMB 
temperature map. 



3 By adding a small e-term to the response W, and thus making it invert- 
ible, the third step also holds generally. 
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Let us now look at 

— {TTcorr^)v(T,P 



rrivec f\ 
uncorr-^corr /V(T,P \ p) 



\P) 



/rpvec rprec f\ 
\ i corr i corr /V(T,P\p) 



= (TP ] ) v{t , p \ p) {Np + WSpW t y 1 WS P ,T 

- S T ,pW^ (Np + WSpW^y^PP^rtp , p) 

(n p + wspW^y V&, T 

S^pW^TVp + WSpW^WSpv 

- St,pW\N p + WSpW^-^Np + ^pVK 1 ") 
(iVp + WSpW^^WSp^ 

St,pW\n p + VF5pVK t )" 1 V^5p ; T 

- St,pW\N p + ^p^)" 1 ^ 

= o, 

where we have inserted eqs dAlK dA2K and dA31 >. QED 



APPENDIX B: PROOF OF VANISHING CORRELATION 
BETWEEN Pcorr AND P u r ^ c corr 

For the splitting of the polarization map, we first prove that the 
unfiltered uncorrected map defined in eq. (|23K P^Zom 1S uncor- 
rected with -Pcorr • 

^uncorr^corr /V(T,P\p) 

— (PPcorr^)v(T,P \ p) ~ W (-Pcorr -Pcorr ^ ) V(T,P \ p) 

— (PT^) V ( T ^ P \ P )S T 1 St,p 

— W Sp^rSrp 1 {TT^)-p(T,p | P )S T 1 St,p 
= W Sp^tSji 1 St,p — W Sp^tSj, 1 StSj, 1 St,p 

— W Sp^tSj, 1 St,p — W Sp^tSj, 1 St,p 

= 0. (Bl) 

From the above, we readily obtain that also the Wiener filtered un- 
corrected map, 

Puncorr = [(Sp-Sp^ST.py'+W^N-'W]- 1 

is uncorrected with P CO rr : 



-)raw 
ancorr 5 



(Puncorr Pcorr )v(T,P \ p) 

[(Sp - Sp^S^St^)' 1 + W^N^W]' 1 

W Np 1 (Puncorr Pcorr }v(T,P\p) 



and independently from each other. We then work out the probabil- 
ity for at least one of the axes of Puncorr being such that the axis of 
the temperature map is included in the la region around it. 

For simplicity, we assume that the the la regions are sym- 
metric circles around the axes, with radius a ~ 45° for both the 
quadrupole and the octopole. The solid angle A spanned by such a 
la region is well approximated by A w 7rcr 2 0The probability of 
at least one of the la regions hitting the axis of evil is just the solid 
angle spanned by the two la regions divided by the solid angle of 
the hemisphere, 2ty. However, the solid angle spanned by the two 
la regions depends on the overlap B between them, it is 2 A — B 
to avoid double counting of the overlapping area. Given the angular 
separation a between the axes of the quadrupole and the octopole, 
the overlap can be computed as follows: 



B(a) 



a arccos 



(CI) 



which can be derived from the geometry of the problem in flat- 
sky approximation. We marginalise the hitting probability over the 
overlap B(a), using the fact th at a is distributed as V{a) — sin(a) 
dde Oliveira-Costa et al.ll2004h : 



P(hit) 



k/2 



/ P(hit | B(a))V(a) da 
B(a) 



2A 



2tt 



(a) da w 50% . 



(C2) 



(B2) 



QED 



APPENDIX C: PROBABILITY FOR CHANCE 
ALIGNMENT IN AN ISOTROPIC UNIVERSE 

We would like to assess whether the rough alignment of the axis 
of the quadrupole in Puncorr actually provides us with some infor- 
mation about the axis of evil. We therefore compute the probability 
for at least one of the axes of Puncorr aligning with the axis of the 
temperature map in an isotropic universe. To this end, let us take 

the preferred axis in the temperature map T as given, and assume 4 This flat-sky approximation differs from the actual value of the solid an- 
that the axes of Puncorr are distributed isotropically over the sky gle by 6 per cent. 
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